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Yogesh N. Joglekar, Clinton Thompson, and Gautam Vemuri 
Department of Physics, Indiana University Purdue University Indianapolis (lUPUI), Indianapolis, Indiana 4^202, USA 

(Dated: January 15, 2013) 

We investigate the single-particle time evolution and two-particle quantum correlations in a 
one-dimensional A^-site lattice with a site-dependent nearest neighbor tunneling function ta{k) = 
to[k{N — k)]^^'^ . Since the bandwidth and the energy levels spacings for such a lattice both depend 
upon a, we show that the observable properties of a wavepacket, such as its spread and the relative 
phases of its const it utents, vary dramatically as a is varied from positive to negative values. We also 
find that the quantum correlations are exquisitely sensitive to the form of the tunneling function. 
Our results suggest that arrays of waveguides with position-dependent evanascent couplings will 
show rich dynamics with no counterpart in present-day, traditional systems. 



Introduction: Idealized lattice models have been popular 
in physics due to their analytical and numerical tractabil- 
ity [1 , the absence of divergences associated with the 
ultraviolet cutoff [U |3], the availability of exact solu- 
tions i4|, and the ability to capture counter-intuitive 
physical phenomena including the bound states in re- 
pulsive potentials [5 . Over the years, these models have 
been successful in describing a diverse array of physical 
systems with bosons, fermions, and quantum spins, with 
short- or long-ranged interactions, such as electronic ma- 
terials, optical lattices ^|7^, and, most recently, evanes- 
cently coupled optical waveguides [8]. The prototypical 
lattice models have a constant nearest neighbor tunneling 
amplitude, and the effects of ubiquitous disorder and im- 
perfections are taken into account via random impurity 
potentials and small, random variations in the tunnel- 
ing amplitude. They are sufficient to capture important 
physical phenomena such as Anderson localization [9]. 

In recent years, coupled optical waveguides have be- 
come a paradigm for the realization of an ideal one- 
dimensional lattice model with tunable tunneling and 
on-site potential, as well as non-Hermitian parity- and 
time-reversal {VT) -symmetric potentials [TOl |TT] . They 
have been used demonstrate several phenomena from 
condensed matter physics and quantum optics, such as 
Bloch oscillations [12], Dirac zitterbewegung [13], Tal- 
bot effect [14 , and quantum random walks [T5^. An- 
derson localization due to random on-site potential, in- 
troduced by a randomly varying refractive index, has 
been experimentally observed in waveguides with a con- 
stant nearest-neighbor tunneling [16 . Two-particle An- 
derson localization, quantum statistics effects, and quan- 
tum and classical correlations have been theoretically ex- 
plored in such waveguides as well [17 . The spontaneous 
T^T-symmetry breaking has been observed in two cou- 
pled waveguides with T^T-symmetric complex index of 
refraction [18]. Most of these cases, with the notable ex- 
ception of Refs. [TQ'.W, have primarily focused on on-site 
disorder effects in a one-dimensional lattice with roughly 
constant nearest-neighbor hopping that, in the contin- 
uum limit, translate into disorder effects on particle with 
a finite mass and a quadratic dispersion; in particular, 
properties of itinerant quantum particles in a lattice with 
position-dependent tunneling amplitude have not been 



extensively explored. 

In this paper, we show that a one-dimensional lat- 
tice with position-dependent tunneling function ta{k) = 
to[k{N — /c)]"/^ = ta{N — k) has a rich dynamics with 
no counterpart in the traditional lattice. This model is 
motivated by the robust T^T-symmetric phase of its non- 
Hermitian counterpart [21 . The results that we report 
here on single-particle propagation and localization, and 
two-particle correlations in such lattices can be investi- 
gated in an array of coupled optical waveguides. Phys- 
ically, one would need to engineer the waveguides such 
that the coupling between adjacent waveguides has the 
form mentioned above. Globally, we find that the pa- 
rameter a has a significant effect on the evolution of the 
input wavepacket, and that tuning a allows one to tailor 
certain aspects of the wavepacket evolution. 

Our salient results are as follows: i) the spread of a 
wavepacket, after propagating a certain distance along 
the waveguides, monotonically depends upon a; in par- 
ticular, when a < — 1, the wavepacket spread is negligible 
for physical propagation distances, ii) when a = 1, the 
phase-information in the initial state of the particle is ac- 
cessible only within windows around certain propagation 
distances, and the size of these windows can be controlled 
by the location of the input waveguide; when a 7^ 1, this 
phase information is, in principle, always accessible, iii) 
for two quantum particles injected into adjacent waveg- 
uides, the two-particle correlation function is exquisitely 
sensitive to a and the location of the input waveguide. 

As we discuss below, these results are not substantially 
affected by a "weak" disorder. They show that coupled 
optical waveguides with specifically engineered tunneling 
functions may provide novel, heretofore unexplored, re- 
alizations of lattice models with tunable energy levels, 
densities of states [21 , and two-particle correlations. 
Tight-binding Model: We consider an array of N waveg- 
uides described by the Hamiltonian with open boundary 
conditions, 

N-l N 

Ha = - ^ ta{i) (4+1^* + ^^Vialai, (1) 

i=i i=i 

where is the creation operator for a particle at site /c, 
ta{k) is tunneling amplitude between sites k and /c + 1, 
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and Vk represents the potential, determined by the local 
index of refraction, at site k. The tunneling amplitude 
t{k) is determined by the evanescent coupling between 
waveguides k and /c + 1, and can be tuned by varying the 
width of the barrier between the two waveguides [19 . 

A Hamiltonian eigenfunction = Xlfe'^fc^Il^) ^^^^ 

energy E'^ satisfies the difference equation 



ta{k-l)i^]:_^^ta{k)^, 



(2) 



where we have considered a constant index of refraction 
riji^ which results in a constant shift in the energy eigen- 
values. The eigenvalue spectrum for Eq. [2] is symmetric 
about zero Hence, the bandwidth of the spectrum, 
defined as the difference between the maximum and min- 
imum eigenvalues, is A^, = £^max - E^in = '^E^a^. 
Note that when a > 0, the tunneling function ta{k) is 
maximum at the center of the waveguide array whereas 
when a < 0, it is maximum at the ends. As a result, 
when N ^ 1 the bandwidth Aa{N) of the Hamilto- 
nian Ha increases monotonically with a. It is natural 
to use the inverse-bandwidth as the characteristic time, 
Ta = 2h/Aa^ and = cTa/riR as the characteristic 
distance along the waveguide where h = h/{27r) is the 
scaled Planck constant and c/ur is the speed of light in 
a waveguide. Note that since Tq, and La are monotoni- 
cally decreasing functions of a, a waveguide array with a 
fixed physical length will correspond to " short-time" sce- 
nario when a < and "long-time" scenario when a > 0. 

Figure [l] shows the time-evolution of a wavepacket that 
is initially localized in waveguide mo = 5 in an array 
of = 25 waveguides. The vertical axis denotes dis- 
tance along the waveguide for a fixed physical length 
of the waveguide L/{hc/nto) = 50. The three vertical 
panels correspond to a = (left), a = 1 (center) and 
a = —1 (right). When a = 0, the traditional model, 
the wavepacket broadens as it travels down the length 
of the waveguide array. When a = 1, the energy levels 
are given by En = ±to(A/' — 1), ±to(A/' — 3), . . .; the level 
spacing is constant and the bandwidth is Aa=i{N) = 
2(A^ - l)to [201 El]. Therefore we obtain perfect re- 
construction of the wavepacket, shown by white cir- 
cles, at mirror-symmetric positions {N + 1 — mo) = 21 
and mo = 5. It should be noted that the wavepacket 
first reconstructs at the mirror symmetric waveguide, i.e. 
21st site, and then alternates between sites 5 and 21. 
For a = —1 (right panel), the wavepacket spread over 
the same distance along the waveguide is significantly 
smaller, consistent with the smaller bandwidth of the 
Hamiltonian. 

To explore the intrinsic a-dependence of the time- 
evolution, in the rest of the paper, we consider waveg- 
uides with the same normalized length L/La = 100; phys- 
ically, this will correspond to waveguides with different 
a-dependent lengths. Figure |2] shows the time-evolution 
of a wavepacket initially at mo = 5 for a = 1 (left panel), 
a = 2 (center panel) and a = —1 (right panel) in an array 
of A" = 25 coupled waveguides. The vertical-axis shows 
distance (time) in the units of La{Ta). Apart from the 
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FIG. 1. (color online) Time evolution in an array with a fixed 
length L/(hc/nRto) — 50 shown along the vertical axis. Top 
three panels show the probability-amplitude time-evolution 
plots of a single photon injected in waveguide mo = 5 in an 
array of A = 25 waveguides with a = (left panel), a — 
(center panel) and a — —1 (right panel). The vertical axis in 
each panel represents the distance along the waveguide. When 
a = the wavepacket, initially localized at mo = 5, spreads as 
it travels along the waveguide. For a = 1, because the energy 
levels are equidistant, the wavepacket is periodically localized 
at mirror symmetric positions (A + 1 — mo) — 21 and mo — 5 
(white circles). When a — —1, the wavepacket spread is 
noticably smaller over the same length of the waveguide. The 
bottom panel shows the dimensionless bandwidth Aa{N)/to 
(vertical) vs. A (horizontal) on a logarithmic scale for 25 < 
A < 250. We see that when a > 1, Ac,(A) oc A°^, whereas 
for a < -1, Ac (A) oc A""^/^. Therefore, a sample with a 
given physical length represents "short-time" evolution when 
a < and "long-time" evolution when a > 0. 



perfect reconstruction at mirror-symmetric points that 
occurs when a = 1, we see that, in contrast to the be- 
havior in Fig. [l] the spread of the wavepacket is qualita- 
tively similar for all a over the normalized length-scales 
(or time-scales). The bottom panel shows, for a = — 1, 
the time-evolution of a single photon injected near the 
edge, mo = 2; the horizontal axis shows the normalized 
distance (time). In this case, the photon remains at the 
edge due to localized edge eigenstates that are generically 
present when a < [2T] . 

Lastly, we consider a wavepacket that is initally local- 
ized in mirror-symmetric waveguides, l^/^^) = (|mo) + 
g«0|7V + 1 — mo))/V2. We obtain the time-evolution 
of the ^-dependent probability amplitude A(0, k) = 
\{k\ip(i){t))\ and use the maximal difference A{t,k) = 
A{0,t,k) — A{7r/2,t,k) as the indicator of the phase in- 
formation. Note that since the initial state is localized 
in two spatially separated regions, information about the 
phase (p will become visible in A{t, k) only after a time 
when the partial waves from the two mirror-symmetric 
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FIG. 2. (color online) Top panels show the time evolution of a 
photon injected at mo = 5 in an array of A/" = 25 waveguides 
with a = 1 (left), a = 2 (center), and a = — 1 (right). The 
wavepacket spread is similar over normalized length-scales. 
The bottom panel shows, for a = —1, the time evolution 
of a photon injected near the edge, mo = 2. The strong 
localization of the photon near the edge is due to the presence 
of localized edge eigenstates that occur when a < 0. 

sites interfere with each other. Figure [s] shows A{t^ k) 
for an array of TV = 25 waveguides with mo = 1, where 
the vertical axis represents normalized distance (time) 
along the waveguide. When a = (left panel) and 
a = —1 (right panel) the phase information, indicated by 
a nonzero value of A, is visible at all times, as is expected 
for a clean system. When a = 1 (center panel), however, 
the phase information is available only in a restricted 
window in the (t, k) space. The size of this window in- 
creases with mo- Thus, when a = 1, the information 
about the initial relative phase remains inaccessible over 
a large fraction of the parameter space. 
Two-particle correlations: We now explore the effects 
of the tunneling function ta{k) on the two-particle 
(number) correlation function defined by r^^{t) = 
{al^{t)al^{t)an{t)ajn{t)). This function encodes the 
Hanbury-Brown-Twiss quantum correlations in coinci- 
dence detections in waveguides m and n [23]. For an 
initial state where the two particles are localized at sites 
(mo, no), the correlation function becomes 

mn 

it) = \G 

mno 

it)G Wl • (3) 

where Gpq{t) = [ex.p{—iHat/h)]pq is the time-evolution 
operator and ± signs correspond to bosons and fermions 
respectively. When a = 0, the traditional model, proper- 
ties of this correlation function and its dependence on the 
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FIG. 3. (color online) Maximal amplitude difference A{t, k) 
for an array of A/" = 25 waveguides with initial state l^^^) = 
(|1) + e''^\N))/V2. In ah cases, ^ = at short distances 
I /La < 15 where the partial waves from initial waveguides 
do not interact with each other. When a = (left panel) 
and a = — 1 (right panel) the phase information persists, as 
is expected for a clean system. When a = 1 (center panel), 
however, the phase information is accessible only in certain 
(diamond shaped) windows. 



initial state have been extensively investigated [17] . Since 
T^^{t) is determined by the time-evolution operator, it 
follows that the bosonic and fermions correlations will be 
qualitatively different when a ^ 0. In particular, when 
a = 1, the constant energy level spacing implies that 
^mri^i'^) is periodic in time or, equivalently, in the dis- 
tance along the waveguide; since the maximum spread of 
a wavepacket initially confined at position 1 < mo < N/2 
is approximately 2mo, it follows that the spatial extent 
and shape of the correlation function in the (m, n) plane 
can be controlled by appropriate initial conditions. 

Figure [4] shows FJ^^ for an array with = 40 waveg- 
uides and (mo, no) = (1,2). The left panels show the 
results for a = 1 for bosons (top) and fermions (bottom) 
at time t/T^ = 25. In contrast to the a = case [17], 
the correlation function is strongly localized at all times, 
and has only two peaks with a single nodal line sepa- 
rating them. The right panels correspond to a = 2 and 
t/To; = 55. At this time, the bosonic correlation func- 
tion (top) is localized near the second edge, with a nearby 
parabolic nodal region. On the other hand, the fermionic 
correlation function (bottom) is sharply localized in one 
direction and extended in the other, with a broad nodal 
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FIG. 4. (color online) Left panels show the correlation func- 
tion r^^(t) for an array of = 40 waveguides at t/T^ = 25 
and a = 1. The initial state of the system has two particles in 
the first two waveguides. The correlation function remains lo- 
calized, and develops only two symmetric peaks with a single 
nodal line. Right panels correspond to a = 2 and t/To, = 55 
with the same initial conditions. We see that the bosonic 
correlation function (top) is localized near the second edge, 
whereas its fermionic counterpart (bottom) is localized in one 
direction and extended in the other direction. 

region around the diagonal. (We recall, from the central 
panel in Fig. |2| that when a = 2, a wavepacket starting 
near the edge localizes substantially near the other edge 
when t/Tc, ~ 50.) These results show that the quantum 
statistics lead to nontrivial correlations for a = 2 case 
that are dramatically different from the a = case [17] 
or the a = 1 case. 

Discussion: In this paper, we have shown that modifying 



the tunneling function in a tight-binding Hamiltonian, 
which can be realized by an array of coupled waveguides, 
produces a wide range of wavepacket evolutions that are 
not seen in traditional models. The tunneling function 
toc{k) affects the wavepacket properties through the band- 
width Aq, and energy level spacings, both of which are 
dependent on a. 

For waveguides with a fixed length, we have shown that 
there are qualitative differences in the wavepacket time- 
evolution depending on whether a is positive or negative. 
For example, when a = 1, the equidistant energy lev- 
els lead to periodic behaviors such as wavepacket recon- 
struction [20 ; when a < 0, a wavepacket near the edge 
remains localized due to edge eigenstates. In addition, 
we have shown that when a = 1 the phase-information 
about an initial state remains inaccessible over a large 
region of the parameter space, whereas when a ^ 1, it is 
accessible. 

We have shown that the tunneling function modifies 
quantum correlations in a non- trivial manner. For ex- 
ample, when a = 1, the size and the shape of bosonic 
and fermionic correlations can be tuned by the choice of 
initial waveguides; the periodicity of these correlations 
follows from the equidistant energy spectrum. For the 
same initial conditions, when a = 2, we find that the cor- 
relations, including the no-coincidence region for bosons 
and fermions, are dramatically different. 

These result are applicable for a clean, disorder-free 
system. For a finite lattice, a weak disorder Vd/ <C 1, 
will localize a wavepacket to its initial waveguide [16] 
after a sufficiently long time Ti ^ Tq,, or distance along 
the waveguide. The disorder and propagation-distance 
thresholds, as well as the effect of a weak nonlinearity, 
however, depend upon a [24]. Our results, thus, remain 
valid at times ^t<^Ti. 
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